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We investigate the dynamical charge-charge correlation function and the optical conductivity in 
weakly doped antiferromagnets using Mori-Zwanzig projection technique. The system is described 
by the two-dimensional t-J model. The arising matrix elements are evaluated within a cumulant 
formalism which was recently applied to investigate magnetic properties of weakly doped antiferro- 
magnets. Within the present approach the ground state consists of non-interacting hole quasiparti- 
cles. Our spectra agree well with numerical results calculated via exact diagonalization techniques. 
The method we employ enables us to explain the features present in the correlation functions. We 
conclude that the charge dynamics at weak doping is governed by transitions between excited states 
of spin-bag quasiparticles. 



The dynamical properties of the normal state of the superconducting cuprates are still not completely understood. 
A much discussed question about strongly correlated electronic systems in two dimensions is whether they can be 
described as Fermi liquids or not. In Fermi liquids the low- lying spin and density excitations can be considered as 
particle- hole excitations of fermionic quasiparticles. In contrast, one-dimensional correlated systems show spin-charge 
separation, i.e. spin and charge excitations are different elementary excitations, and the "physical" electron can be 
interpreted as a composite object. 

Theoretical progress in this field is mostly based on numerical techniques, especially on exact diagonalization 
methods and Lanczos calculations However, the system sizes presently accessible leave many problems 

unresolved. Meanwhile, the spin response of doped antiferromagnets has been studied in a number of analytical 
papers but only few authors have investigated the charge density response function . 

In this paper we present analytical results based on projection technique calculations for dynamical charge-charge 
correlation functions in weakly doped 2D antiferromagnets. The system is described by the 2D t-J model: 

H = -tY. + clh.) + JY1 (S.S, - ^) . (1) 

The symbol (ij) refers to a summation over pairs of nearest neighbors. The two antiferromagnetic sublattices will be 
denoted by t and |. The electronic creation operators exclude double occupancies. 

We study the Laplace transform of the time- and wave- vector dependent charge-charge correlation function at zero 
temperature defined by 

G,{k,z) = {^PoIpI^pM (2) 

Z Ij 

Here, pk = cr "^k+q (T<^qcr — J2icr ^iij^ia ^hc Fouricr-transformcd charge density operator. |?/'o) denotes the 
exact ground state of the system, z = a; + 177 is the complex frequency variable. The Liouville operator L is a 
superoperator defined by LA = [H,A]^ for any operator A. At zero temperature the optical conductivity cr(cij) for 
frequencies a; > is related to the charge response function Gp(k,uj) as follows: 

Re.H = -hm ^ (3) 

where N is the total particle number. Eq. (^) follows from the continuity relation Lpk = k- jk where jk is the charge 
current operator. 

In the following, we apply a projection technique |Q| approach to determine the charge-density response function. 
The arising matrix elements are evaluated using a cumulant formalism which has been recently introduced 
to investigate ground-state properties of correlated electronic systems. We are interested in calculating dynamical 
correlation functions for a set of operators 

G,p{z) = i^olSBl^SBJ^Ijo) (4) 
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with SB^ = B,y — (■0o|-Bi/|'0o)- Using cumulants these correlation functions can be rewritten as |^ 

G,^{z) ^ (cl^oin^ Bl(^-^B^^ n\<j>or (5) 

The operator has similarity to the so-called wave operator (or Moeller operator known from scattering theory). 
Within cumulants it transforms the ground state |0o) of the unperturbed system Hq into the exact ground state \tPq) 
of H = Ho + Hi. Explicitly it is given by [|| 

x^o X - [Lo + Hi) 

The brackets {(f>o \ ... \(t>o)'' denote cumulant expectation values with \(f)o). The dot • in Eq. (^) indicates that the 
quantity inside (...)' has to be treated as a single entity in the cumulant formation. Lq is the Liouville operator 
corresponding to Hq, i.e., LqA = [Ho,A]-. The relation (||) can be applied to either weakly or strongly correlated 
systems because its use is independent of the operator statistics , i.e., it is valid for fermions, bosons or spins. 

Using Mori-Zwanzig projection technique Q one can derive a set of equations of motion for the dynamical correlation 
functions G,y^(z). Neglecting the self-energy terms as explained below it reads: 

Xr\v and LOrfv are the static correlation functions and frequency terms, respectively. They are given by the following 
cumulant expressions: 

^r,u = Y.^U^^BI(LB;)-^\^oYxII. (8) 

A 

X^/i is the inverse matrix of x^ii- These terms describe all dynamic processes within the subspace of the Liuoville 
space spanned by the operators By. 

Next we outline the description of the ground state of the t-J model at weak doping within the cumulant formalism, 
details have been published recently (5 denotes the hole concentration away from half filling, i.e., the system with 
iV lattice sites possesses M = 5iV dopant holes. The Hamiltonian is decomposed into Hq and H\ as follows: 

Hq = Hjsing , Hi = Ht + H± . (9) 

The ground state |(/>o) of the unperturbed Hamiltonian Hq is an antiferromagnetically ordered Neel state with M 
holes. The holes have fixed momenta k„i and are located on the sublattice cr„i {cFm ^Ti i) 

0o} = CkicTi . . . CkMO-A/ Ift'Neel} , 

Within the cumulant method we employ an exponential ansatz [0 for the wave operator fi, i.e., fl = . The 
operator S describes the effect of the perturbation Hi onto the unperturbed ground state {(pa). We are basically 
interested in the charge dynamics of the system. Therefore we assume that background spin fluctuations can be 
neglected. To treat hole motion processes induced by Ht the concept of path operators [|2| is used which leads to the 
picture of spin-bag quasiparticlcs ]l3t . Here we define path concatenation operators A„_p(j) where i denotes a lattice 
site, n the path length, and p the individual path shape: A„^p(i) operating on a hole at site i in the state (^0|) moves 
the hole n steps away and creates a path or string of n spin defects attached to the transferred hole. Note that there 
is a number to„ of different path shapes for a given path length n. For n = 1 there are mi = 4 different paths, for 
n = 2 7712 = 12 and so on. Having defined the excitation operators A„_p(i), the wave operator Q of the cumulant 
formalism takes the form 

( "ma J m„ \ 
^ ^K,pAn,p I , An^p = ^An^p{i) (11) 
n=l p=l J i 
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with parameters \n,p yet unknown. Note the additional definiton Af){i) = i Cjg.cjg. which is the only "path" with 
zero length (p — 1). The path operators ^n,p commute with each other because they only contain spin-flip operators 
destroying Neel order. 

Following ref. ||l^ one attains a non- linear set of equations for the ground-state energy and the coefficients A„.p: 

Eo = (0o|i^^^|0o)^ = {cf>o\AlpHn\cj,aY (12) 



with Q. given by (11). The cumulant expectation values have to be taken with respect to the unperturbed ground 
state (^0|). The set of equations (^2|) together with the additional approximation of independent hole quasiparticles 
leads to a generalized eigenvalue problem. For details see ||lO|. The used picture of independent quasiparticles is 
appropriate for small hole concentrations, i.e., for a dilute gas of holes moving in an antiferromagnetic background. 

For the application of Mori-Zwanzig projection technique one has to choose a set of relevant operators Bv Here 
we are going to neglect the self-energy terms obtained within projection technique, cf. Therefore we need to use 
a sufficiently large set of operators to cover the charge dynamics. One variable to include is the charge density 
operator itself. The other variables should provide a coupling between the ground state ri|0o} a-nd all excited states 
of H. These variables are obtained by applying the perturbation Hi to pk- Assuming that the charge dynamics is 
mainly carried by spin-bag quasiparticles we can define 

i3„,,(k) = ^e'''^' A„,p(z)Ao(i)t = ^e'''^' i?„,p(z) . (13) 

i i 

As explained above, A„^p(i) is a path operator of length n with the path shape p acting on a hole at site i. The 
indices n, p replace the general index v of the dynamical variables B^, in Eqs. (^,||). The Fourier-transformed operator 
Bn.p{i) couples to a hole (at site i) and adds the path A„ p. So the relevant part of the Liouville space spanned by 
the operators -Bn,p contains all individual path states connected to the holes. 
The first of these operators is the hole density operator itself: 

Bo,i(k) - ip(j°'^ = ^E^''^'''^-^ (14) 

The quantity we are interested in is therefore the diagonal correlation function Gi^^O^, z) with i> ~ (0, 1). Note that 
particle density and hole density are equivalent quantities. 
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FIG. 1. Left panel: charge-charge correlation function obtained from the present calculation for t/J = 2.5, very small 5, 
rimax = 5, and different momentum transfers, plotted with Lorentzians using an artificial linewidth of 0.1 t. Right panel: Exact 
diagonalization data from ||l] for comparison, t/J = 2.5 (see also [^). The numerical data have been calculated for a 4 x 4 
periodic cluster and a hole concentration of (5 = 25%. 
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Our results for the charge-response function ImG'p(k,w) at t/J = 2.5 are shown in Fig. 1. The intensity of the 
spectra is proportional to the hole concentration 5. For small momentum the spectral weight is mainly concentrated 
in a peak near lo — Q (with an energy scaling with v^? ■ k for k — > 0). With increasing momentum we find a transfer 
of spectral weight to higher energies. Calculations for other values oi t/J show that the broad structures observed 
in the spectra for large momenta scale with t, e.g., the maximum spectral weight in the charge-response function at 
(tt, tt) remains at energies of about 6t. The right panel of Fig. 1. shows exact diagonalization data for S — 25% taken 
from ref . |^ . We observe a good agreement of the spectra, some differences may be either due to the small number 
of sites in the numerical calculations (which leads to an energy gap in all spectra and other finite-size effects) or due 
to the neglection of relaxation processes in the present calculation (such processes would produce finite linewidths). 

For the interpretation of the calculated spectra one can consider the one-hole spectrum of the Hamiltonian within our 
approximation. From the diagonalization of the one-hole problem in the subspace of path operators A^p one obtains 
several bands for the spin-bag quasiparticles. As is well known, the lowest band has a minimum at (±7r/2, ±7r/2) and 
a bandwidth of about 2 J. It corresponds to a quasiparticle with s-like symmetry. The spin-bag states in the higher 
bands have nodes in the coefficients A„.p. The low-energy peak in the spectra (visible at momentum transfer {tt/2, 0) 
) corresponds to excitations within the lowest quasiparticle band whereas the high-energy part arises from transitions 
to higher bands. 
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FIG. 2. Optical conductivity for small hole concentration and different t/J. The spectra have been calculated from the 
charge-charge correlation function using (^). 



Using Eq. (^) one can deduce the optical conductivity from the calculated charge response function. Results 
for Ke(j{uj) are displayed in Fig. 2. They show again good agreement with the numerical results from ref. ^ for 
6 = 12.5%. The main peak (at low energy) scales with J, it is located at « 1.7 J. It arises from the transition between 
the first and second quasiparticle band, i.e., from the s-like groundstate to a p-like state. So our calculation supports 
the discussion given in ref. j^] where only these two states have been considered in a simplified analytical calculation 
for the main peak. We want to emphasize that taking into account all string states (up to a truncation length) as done 
in the present work reproduces not only this main peak but also the incoherent continuum found at energies up to 6t. 
The low-energy peak corresponding to uj ^ 0.3 eV for J — 0.15-0.2 eV which was also found in numerical studies of 
the t-J and Hubbard models is supposed to coincide with the mid- infrared peak observed in optical spectra of high- Tc 
superconductors (see e.g. Q). Within the present calculation we do not obtain a Drude peak DS{uj) since Eq. (||) 
is valid for non-zero frequency only, i.e., it does not include the diamagnetic part of the current. The different scaling 
behaviour of Gp(k, w) (for large k) and cr(w) with t/ J arises from the internal structure of the spin bag forming the 
quasiparticle. Details will be published elsewhere |15| . 

The agreement of our data with the exact diagonalization results jl]j^ is remarkable because these numerical 
calculations are done for rather large hole concentrations (e.g., 25%) where only short-range magnetic order is present 
in the high-Tc materials. In contrast, the present calculations are based on magnetic long-range order. Thus we 
conclude that whether the quasiparticles move in a long-range ordered background or not does not have an essential 
influence on the charge response of the system. The dynamics is mainly determined by the local antiferromagnetic 
order in the vicinity of the hole quasiparticles i.e., the magnetic correlation length has to be of the order of the 
quasiparticle size. 

One should note that a recent slave-boson approach Q to the charge dynamics in the t-J model also reproduces 
some key features of the charge response function, namely a peak at small energy for small k and a broad featureless 
continuum at energies of several t. The theory presented there is based on a Fermi liquid picture and completely 
neglects antiferromagnetic correlations in the ground state. However, in the high-Tc materials short-range correlations 
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arc present also at higher hole concentrations. We believe them to be important for the structure of the spectra. 

Summarizing, we conclude that the charge dynamics in the weakly doped t-J model can be well explained within 
the picture of independent spin-bag quasiparticles and their excitations. The present calculation may serve as a basis 
for the investigation of the charge dynamics in more realistic models for the cuprate superconductors, as e.g. the 
three-band Hubbard model. 
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